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ABSTRACT 

Circumstellar dust disks have been observed around many nearby stars. However, 
many stars are part of binary or multiple stellar systems. A natural question arises 
regarding the presence and properties of such disks in systems with more than one 
star. To address this, we consider a sample of 449 systems (spectral types A-M) 
observed with the Herschel Space Observatory as part of the DEBRIS program. We 
have examined the stellar multiplicity of this sample by gathering information from 
the literature and performing an adaptive optics imaging survey at Lick Observatory. 
Five new companions were revealed with our program. In total, we identify 188 (42%) 
binary or multiple star systems. The multiplicity of the sample is examined with 
regards to the detection of circumstellar disks for stars of spectral types AFGK. 

In general, disks are less commonly detected around binaries than single stars, 
though the disk frequency is comparable among A stars regardless of multiplicity. 
However, this sample reveals the period distribution of disk-bearing binaries is consis¬ 
tent with that of non-disk binaries and with comparison field samples. We find that the 
properties of disks in binary systems are not statistically different from those around 
single stars. Although the frequency of disk-bearing FGK binaries may be lower than 
in single star systems, the processes behind disk formation and the characteristics of 
these disks are comparable among both populations. 

Key words: binaries: general — infrared: stars — planetary systems: formation. 


1 INTRODUCTION 

Disks rich in gas and dust around young stars are sites of 
planet formation. Jovian-class planets need to form rapidly, 
as gas within the disks dissipates over a period of a few mil¬ 
lion years (Zuckerman et al. 1995; Pascucci et al. 2006). The 
dust in the system is also removed eventually, either through 
accretion onto larger objects, stellar winds, or radiative pro¬ 
cesses (Haisch et al. 2001; Uzpen et al. 2009). However, col¬ 
lisions between any planetesimals in the disk can generate a 
second population of dust. These second-generation systems 
are known as debris disks and generally contain very little 
gas (Zuckerman 2001; Wyatt 2008, and references therein). 
Our own solar system’s Kuiper belt may be analogous to 
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these debris disks (e.g., Luu & Jewitt 2002, and references 
therein). The study of disks, from the gas-rich protoplan¬ 
etary disks to the gas-poor debris disks, is necessary for a 
more complete understanding of the formation and evolu¬ 
tion of planetary systems. 

Approximately half of all stars are in binary or mul¬ 
tiple star systems (Duquennoy & Mayor 1991; Eggleton & 
Tokovinin 2008; Raghavan et al. 2010; Duchene & Kraus 
2013). Given how common binary stars are, it is important 
to address the properties of planet formation in such com¬ 
mon systems. About 20% of known extra solar planets re¬ 
side in wide binaries with separations of order 100s of AU 
(Raghavan et al. 2006; Eggenberger et al. 2007). The Ke¬ 
pler satellite has also revealed several planets orbiting pairs 
of close binary systems (e.g., Doyle et al. 2011; Welsh et 
al. 2012; Orosz et al. 2012a,b; Schwamb et al. 2013; Kos- 
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tov et al. 2014). Earlier work with eclipse timing variations 
have suggested planets around the binaries HW Virginis, 
CM Draconis, and NN Serpentis (a post-common envelope 
binary) (Deeg et al. 2008; Lee et al. 2009; Beuermann et al. 
2010 ). 

Circumstellar disks can be used as indirect evidence of 
planets in binary systems. For example, the close (3.4-day) 
main sequence binary BD +20 307 displays a large amount 
of warm dust in the terrestrial planet zone (Song et al. 2005), 
which can be interpreted as the result of a planetary-scale 
collision in this ~1 Gyr old binary system (Zuckerman et al. 
2008; Weinberger et al. 2011). Disk studies among young 
pre-main sequence systems have shown that the presence 
of a nearby companion star can readily truncate and dis¬ 
perse disks (Jensen et al. 1996; Bouwman et al. 2006; Cieza 
et al. 2009; Andrews et al. 2010; Kraus et al. 2012), in good 
agreement with what is expected from numerical simula¬ 
tions (Lubow & Artymowicz 2000; Artymowicz & Lubow 
1994). Older debris disk binary systems also show a simi¬ 
lar behavior (Trilling et al. 2007; Rodriguez & Zuckerman 
2012). Studying circumstellar and circumbinary disks is then 
a complementary way to explore the properties of planet for¬ 
mation in binary star systems. 

The Herschel Space Observatory 1 (Pilbratt et al. 
2010) offers the best opportunity to date of performing 
a large scale volume-limited survey for debris disks, the 
ideal approach to systematically analyze the interplay be¬ 
tween debris disks and stellar multiplicity. The DEBRIS 
(Disc Emission via a Bias-free Reconnaissance in the In¬ 
frared/Submillimetre; Matthews et al. 2010) program has 
explored nearby star systems to search for infrared excesses 
indicative of disks at wavelengths longward of 70/xm. Prior 
work on the multiplicity of debris disk stars has been lim¬ 
ited by small samples with strong biases; for example, con¬ 
sidering only debris disk systems (Rodriguez & Zuckerman 
2012) or only binary and multiple star systems (Trilling et 
al. 2007). The DEBRIS sample has observed stars regard¬ 
less of prior known disks or stellar multiplicity and thus 
provides a better sample in which to explore the relation¬ 
ships between these two phenomena. Several DEBRIS bina¬ 
ries/multiples with disks have been individually explored in 
detail (Kennedy et al. 2012a,b, 2014). In this paper, we ex¬ 
plore the multiplicity statistics of the DEBRIS sample and 
its role in the detection of circumstellar and circumbinary 
disks. 


2 MULTIPLICITY IN THE DEBRIS SAMPLE 

The DEBRIS sample consists of 451 stars, of which 2 are 
not primaries (Fomalhaut B &C; see Mamajek et al. 2013; 
Kennedy et al. 2014), and includes the 5 systems observed 
by the Herschel Guaranteed Time (GT) disks program (PI: 
Olofsson). The majority of the stars were observed by Her¬ 
schel as part of the DEBRIS program, with some observa¬ 
tions coming from the DUNES program (Eiroa et al. 2010). 
The sample of 449 primaries is roughly equally divided in 
spectral types A through M and is volume-limited for each 

1 Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA. 
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Figure 1 . Distribution of calculated and measured orbital peri¬ 
ods for multiples (gray histogram). Triples and higher order mul¬ 
tiples have more than one period per system. The dashed line is 
the period distribution from Eggleton & Tokovinin (2008) and the 
solid curve is from Raghavan et al. (2010). Both are normalized 
to the number of periods in the DEBRIS sample. The blue shaded 
histogram shows the distribution of periods for disk-bearing mul¬ 
tiples in our sample. 


individual spectral type, where the limit is 45.5, 23.6, 21.3, 
15.6, and 8.6 pc for A, F, G, K, and M spectral types, re¬ 
spectively (Phillips et al. 2010). Herschel observations with 
PACS (Poglitsch et al. 2010) and SPIRE (Griffin et al. 2010) 
have been performed for the DEBRIS sample in order to 
search for and characterize far-IR emission from circum¬ 
stellar (or circumbinary) dust. Except for cases in which a 
companion fell within the field of view, we didn’t explicitly 
aim to observe companions with exception of Fomalhaut’s 
companions, of which companionship was unknown until re¬ 
cently (see Mamajek et al. 2013; Kennedy et al. 2014). Thus, 
we may be missing some circumsecondary disks and thus 
focus on circumprimary or circumbinary disks only. Simul¬ 
taneously to these Herschel observations, we have gathered 
literature data and adaptive optics (AO) observations (see 
Section 3) in order to characterize the multiplicity of stars 
in this sample. 

Given that the majority of the stars in the DEBRIS 
sample are close to Earth and well-studied in the literature, 
published data concerning stellar multiplicity exists for a 
large fraction of the sample. Our AO observations comple¬ 
ment this literature search (see Section 3). We present the 
full list of stars and multiplicity information in Table 1. Some 
systems have tentative evidence of being astrometric bina¬ 
ries, but no additional information exists to confirm this. 
These systems are treated as single stars for this study. Bi¬ 
naries and higher multiples are listed with more details in 
Table 2. In general, we only have measurements of the pe¬ 
riod of spectroscopic binaries and projected separation for 
visual binaries. We derive the remaining quantities (ie, pe¬ 
riod or semimajor axis) assuming the measured separation is 
the orbital semimajor axis, orbits are circular, and adopting 
main-sequence masses from Baraffe et al. (1998) and Siess et 
al. (2000). Which quantity is derived is indicated in Table 2. 
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Table 1. DEBRIS Sample 


Name 

UNS ID“ 

Multiple? 

Notes 

References 

HD 38 

K050 

Y 

(8.2K6+9.9M2; 6.041”) 

HIP, SHAH 

HD 166 

G030 

N 

optical double; LAF07: background source 
10.23” away 

WDS, LAF07, KIY08 

HD 693 

F069 

N 

F7V 

ET08, TAN10 

HD 739 

F096 

N 

F4V 

ET08 

HD 1237 

G070 

Y 

(G6+M4; 3.857”) 

CHA06 

HD 1326 

M011 

Y 

(8.31M2V+11.36M6V; 2600y, 41.15”) WDS: 
AC pair ruled out; TAN10: companions to A 
& B found, but no common motion 

WDS, VB6, TAN 10 

HD 1404 

A103 

Y 

A2V, proper motion companion found in Lick 
AO images 

ET08, ThisWork 

HD 1581 

F005 

N 

F9V 

ET08 

HD 1835 

G118 

N 

optical double, different proper motions 

WDS, RAG 10 

HD 2262 

A021 

N 

A7V 

ET08 

HD 3196 

F092 

Y 

(5.61(F8V + ?; 2.082d) + 6.90G0V; 6.89yr 
e=0.76) 

ET08 

HD 3443 

G044 

Y 

(6.37G8V + 6.57G9V; 25.09yr e=0.24) 

ET08 

HD 3651 

K045 

Y 

(4.55+16.87T7.5; 42.9”) WDS: AB optical, AC 
common proper motion 

WDS, KIY08, LUH07 

HD 4391 

G041 

Y 

(5.8G1V + (13.5; 16.6 251AU) + (14M5; 
49 742AU)); ET08 incorrectly lists as single, 
WDS entries somewhat contradictory 

WDS, RAG 10 

HD 4628 

K016 

N 

K2V 

ET08 

HD 4676 

F124 

Y 

((F8V + F8V; 13.82d e=0.24) + ?; 0.25”) 

ET08 

HD 4747 

G089 

Y 

(?+?; 6832d e=0.64) 

SB9 

HD 4813 

F038 

N 

F7IV-V 

ET08, TAN10 

HD 4967 

K127 

Y 

(8.16K5+15.2M5.5; 16.8”) 

WDS, POV94, 

HAW97, GOU04 

HD 5133 

K089 

N 

no stellar companion 

RAG 10 

HD 5448 

A090 

N 

A5V 

ET08 


Systems in the DEBRIS sample. The notation follows that of Eggleton Sz Tokovinin (2008). 
The full table is available online. Details on binary and multiple systems are listed in 
Table 2. Systems with multiplicity flag “A” have some indication that they are astrometric 
binaries, yet no other information is available. 
a : Unbiased Nearby Stars from Phillips et al. (2010) 

Table references as follows: HIP : Hipparcos Catalog (ESA 1997), SB9 : The ninth catalogue 
of spectroscopic binary orbits (Pourbaix et al. 2004), VB6 : Sixth Catalog of Orbits of Visual 
Binary Stars (Hartkopf et al. 2001), WDS : Washington Double Star Catalog (Mason et al. 
2001), BAI10 : Baines et al. (2010), BAR09 : Bartlett et al. (2009), BERIO : Bergfors et al. 
(2010), CAB07 : Caballero (2007), CAR05 : Carson et al. (2005), CAR11 : Carson et al. 
(2011), CHA06 : Chauvin et al. (2006), CLO07 : Close et al. (2007), DEL99 : Delfosse et al. 
(1999), DER11 : De Rosa et al. (2011), DER12 : De Rosa et al. (2012), DM91 : Duquennoy 
Sz Mayor (1991), EGG07 : Eggenberger et al. (2007), EHR10 : Ehrenreich et al. (2010), 
EIS07 : Eisenbeiss et al. (2007), EKE08 : Eker et al. (2008), ET08 : Eggleton Sz Tokovinin 
(2008), FRA07 : Frankowski et al. (2007), FUH14 : Fuhrmann et al. (2014), GOL06 : Goldin 
Sz Makarov (2006), GOL07 : Goldin Sz Makarov (2007), GOU04 : Gould Sz Chaname (2004), 
HAR12 : Hartkopf et al. (2012), HAW97 : Reid et al. (1997), HIN02 : Hinz et al. (2002), 
HIN10 : Hinkley et al. (2010), JAO03 : Jao et al. (2003), KIY08 : Kiyaeva et al. (2008), 
KOH12 : Kohler et al. (2012), KON10 : Konopacky et al. (2010), LAF07 : Lafreniere et 
al. (2007), LAG06 : Lagrange et al. (2006), LAW08 : Law et al. (2008), LEC10 : Leconte 
et al. (2010), LEI01 : Leinert et al. (2001), LEI97 : Leinert et al. (1997), LEP07 : Lepine 
Sz Bongiorno (2007), LUH07 : Luhman et al. (2007), MAK05 : Makarov Sz Kaplan (2005), 
MAM10 : Mamajek et al. (2010), MAM13 : Mamajek et al. (2013), MAR14 : Marion et al. 
(2014), MAS98 : Mason et al. (1998), MCA93 : McAlister et al. (1993), MCC04 : McCarthy 
Sz Zuckerman (2004), MET09 : Metchev & Hillenbrand (2009), NAK94 : Nakajima et al. 
(1994), NAK95 : Nakajima et al. (1995), NID02 : Nidever et al. (2002), NIE10 : Nielsen Sz 
Close (2010), OPP01 : Oppenheimer et al. (2001), PHI10 : Phillips et al. (2010), POT02 : 
Potter et al. (2002), POV94 : Poveda et al. (1994), RAG10 : Raghavan et al. (2010), ROB11 
: Roberts et al. (2011), SCH07 : Schroder & Schmitt (2007), SCH11 : Schneider et al. (2011), 
SEG00 : Segransan et al. (2000), SHAH : Shaya Sz Oiling (2011), TAN10 : Tanner et al. 
(2010), TOK08 : Tokovinin (2008), TOK10 : Tokovinin et al. (2010), TOK11 : Tokovinin 
(2011), TOK92 : Tokovinin (1992), WIL01 : Wilson et al. (2001). 
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Table 2. DEBRIS Multiples 


Name 

Component 

Period (d) 

Sep (AU) 

Eccentricity 

Derived 0 

Dust R (AU) 

HD 38 

AB 

1.95E+05 

68.2 


p 


HD 1237 

AB 

1.91E+05 

67.5 


p 


HD 1326 

AB 

9.50E+05 

147.3 


N 


HD 1404 

AB 

1.10E+06 

272 


P 

26.4 

HD 3196 

AB 

2.082 

0.035 


a 



AB-C 

2.52E+03 

3.95 

0.76 

a 


HD 3443 

AB 

9.16E+03 

8.93 

0.24 

a 


HD 3651 

AB 

3.72E+06 

474 


P 


HD 4391 

AB 

1.23E+06 

252 


P 



AC 

6.26E+06 

743 


P 


HD 4676 

AB 

13.82 

0.126 

0.24 

a 



AB-C 

4.40E+03 

5.86 


P 


HD 4747 

AB 

6.83E+03 

7.23 

0.64 

a 


HD 4967 

AB 

1.72E+06 

262 


P 


HD 7439 

AB 

9.13E+06 

1162 


P 



Measurements of components of binaries and multiples in the DEBRIS sample. The location 
of blackbody dust grains for dusty systems is also listed. For HD 216956, dust R is listed 
first for A, then C. For HD 223352, the AC pair’s dust is located around the C component. 
The full table is available online. For clarity, we have rounded periods to 1 day when larger 
than 100 days. Similarly, separations in AU are rounded to 0.1 AU when larger than 20AU 
and to 1 AU when larger than 200 AU. We refer readers to the references listed in Table 1 
for the most exact values. 

Quantity derived: P, period; a, semimajor axis; N, none. See Section 2 for details. 


Table 3. Multiplicity Fractions 


Sp. Type 

Number 

This Work 

From RZ12 

From ET08 

Others 

B 



83tL 



A 

35/86 

41 ± 5 

26 ±1 

46.0 

>50 

F 

50/94 

53 ±5 

12 ±1 

47.4 


G 

43/90 

48 ±5 

00 + 15 
°°-10 

45.0 

44 

K 

37/91 

41 ±5 

40+ 21 

^ u -16 

29.1 


M 

23/88 

26±® 



26 

All 

188/449 

42 ±4 

25 ±4 

42.8 



Fraction of multiple stars (all as percentages) broken down by spectral type of the primary 
star with comparisons from the literature (RZ12: Rodriguez & Zuckerman 2012; ET08: Eggle- 
ton & Tokovinin 2008; Raghavan et al. 2010, Duchene & Kraus 2013). Note that the Rodriguez 
& Zuckerman (2012) sample only contains disk-bearing systems. Errors are binomial errors 
estimated as described by Burgasser et al. (2003). 


While clearly not all orbits are circular and projection effects 
have not been incorporated in our analysis (corrections for 
these are small, of order ~10%; see Dupuy & Liu 2011), 
this nevertheless provides useful information in a statistical 
manner. 

The sample has 188 (42%) star systems where two or 
more stars are present. As previously mentioned, we calcu¬ 
late or adopt literature values for parameters of the system, 
such as the orbital period and semi-major axis. In Figure 1, 
we show the period distribution for all multiples in the sam¬ 
ple, including those determined from our AO observations 
(Section 3). This includes all periods, so, for example, triples 
have two periods counted. Also shown is the period distribu¬ 
tions from Raghavan et al. (2010) and Eggleton & Tokovinin 
(2008) normalized to the same number of periods. These 
two latter distributions sample solar-type stars and bright 
systems and are representative of what we would expect for 
our sample. The distributions are remarkably similar to each 


other suggesting that our multiplicity survey is not biased 
against any particular range of periods. 

The multiplicity fraction among A-K stars in our sam¬ 
ple is 40-50%, consistent with some prior studies (Eggleton 
& Tokovinin 2008; Duquennoy & Mayor 1991). However, 
Duchene & Kraus (2013) estimate a lower limit to the mul¬ 
tiple star fraction among A stars of 50% suggesting we may 
still be missing some A-star binaries. Among M-dwarfs we 
find a multiplicity fraction of 26%, in agreement with the 
multiplicity fraction listed in Duchene & Kraus (2013). How¬ 
ever, we note that these faint stars have not been as intensely 
studied for binarity as more massive stars in the DEBRIS 
sample. Furthermore, the DEBRIS Herschel sensitivity to¬ 
wards disks around M-dwarfs is also low (see Figure 2 in 
Matthews et al. 2014). The combination of low detection 
rates (multiplicity or debris disks) results in small number 
statistics whose robustness is too limited for a detailed anal¬ 
ysis. Therefore, we primarily focus on the A-K sample in this 
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article. We summarize the multiplicity statistics, divided by 
spectral type, in Table 3. 


3 ADAPTIVE OPTICS SEARCH FOR 

COMPANIONS 

3.1 Survey Setup and Observations 

We have utilized the adaptive optics (AO) camera IRCAL 
(Lloyd et al. 2000; Perrin et al. 2008) at Lick Observatory 
to search for companions separated from the primary by 
~10-1000 AU. The NIR camera offers a ~0'.'077 pixel scale 
(which provides Nyquist sampling for diffraction-limited ob¬ 
servations at 2.2 micron on the Shane 3m telescope) and a 
20" FOV. Precise pixel scale (including the known slight 
anamorphism of the camera) and orientation was deter¬ 
mined by observations of multiple known binaries. The pre¬ 
cision is estimated to be 1% for the pixel scale and 0.7deg for 
the absolute orientation. The dither pattern typically used 
was a 5"-on-the-side square, giving us full coverage out to 
12'/5 from the primary, and we only consider companions 
within this radius in the survey. Tables 4 and 5 list our ob¬ 
servations and measured parameters for binary and single 
stars, respectively. 

We carried out our observations on various nights in 
2009 June & October, 2010 August, and 2012 March & 
September. A total of 221 DEBRIS targets were observed 
over this time period. This corresponds to 75% of all targets 
with declinations larger than -10 degrees. For our observa¬ 
tions, we used a dithered sequence to remove artifacts and 
cosmic rays and observed with either the Ks or Br -7 filters. 
Objects that had a close companion or appeared extended 
were subsequently observed at J and either H or Fell in or¬ 
der to estimate colors and spectral types. The choice of Fell 
and Br -7 over H and Ks was due to saturation on bright 
targets. The Fell and Br -7 filters have central wavelengths 
of 1.644 and 2.167A with bandwidth FWHM of 0.016 and 
0.020A, respectively. We assume the flux ratio between the 
primary and any companion at Br -7 is comparable to that in 
Ks (and similarly for Fell and H). To confirm companions we 
either obtained a second epoch and tested for common mo¬ 
tion, or estimated the spectral type of the companion with 
our color information and verified that the photometric dis¬ 
tance agrees with that of the primary. We achieve a typical 
contrast of 6 magnitudes (5-cr) at separations >1 arcseconds 
(see Figure 2). The Siess et al. (2000) models, for ages up 
to 1 Gyr, predict a flux ratio of 9-9.5 mag between A0 stars 
(2.5-3 Mq , depending on age) and O.IMq stars at K band. 
Our median detection limit beyond ~2" is 8 magnitudes sug¬ 
gesting we are close to being complete for stellar mass com¬ 
panions. While a few lower mass companions may remain 
to be found, particularly around the later spectral types, 
we do not anticipate that these would significantly change 
the results. Our AO results are incorporated in Figure 1 in 
which we see no strong bias or incompleteness against any 
particular range of orbital periods. 

3.2 Newly Identified Companions 

Most targets observed either had no detected companion or 
were binaries previously known in the literature. However, 


Table 4. IRCAL Companion Measurements 


Name 

Filter 

Sep 

(") 

PA 

(deg) 

Am 

(mag) 

UT Date 

Depth 

(mag) 

HD 38 

Ks 

5.62 

5.87 

0.071 

2009-10-28 

8.47 

HD 1404 

Ks 

6.60 

144.40 

6.682 

2009-10-28 

8.59 


BrG 

6.56 

144.23 

6.746 

2012-09-28 


HD 3196 

Ks 

0.26 

248.43 

0.770 

2010-08-03 

7.97 

HD 16160 

BrG 

1.73 

318.49 

5.142 

2010-08-05 

9.13 

HD 16765 

Ks 

3.82 

300.81 

1.880 

2010-08-04 

9.06 

HD 16970 

Ks 

2.21 

297.83 

1.312 

2010-08-04 

8.93 

HD 19994 

BrG 

2.12 

203.15 

2.830 

2009-10-28 

8.36 

HD 56537 

BrG 

9.29 

35.72 

3.953 

2009-10-28 

8.42 

HD 56986 

BrG 

5.35 

228.90 

3.117 

2009-10-28 

8.7 

HD 76943 

BrG 

0.50 

233.34 

1.161 

2012-03-10 

9.34 

HD 78154 

BrG 

4.02 

348.12 

2.337 

2012-03-10 

9.63 

HD 82328 

BrG 

2.59 

145.18 

5.777 

2012-03-10 

10.37 

HD 82885 

BrG 

6.65 

61.00 

3.753 

2012-03-10 

9.03 

HD 98231 

BrG 

1.50 

196.14 

0.271 

2012-03-10 

9.48 

HD 100180 

Ks 

14.66 

329.39 

1.439 

2012-03-10 

10.32 

HD 101177 

Ks 

8.79 

248.08 

0.879 

2012-03-10 

10.12 


A sample of measurements of companions detected in the 
IRCAL FOV. The depth is the 5 -a limit reached at 4-5" range, 
where the contrast is highest. The full table is available online. 


Table 5. Single Stars Observed with IRCAL 


Name 

Filter 

UT Date 

Depth 

(mag) 

HD 166 

BrG 

2009-10-28 

9.12 

HD 693 

Ks 

2010-08-03 

8.29 

HD 1326 

Ks 

2009-10-28 

8.75 

HD 1835 

Ks 

2010-08-05 

9.07 

HD 3651 

BrG 

2009-10-28 

7.96 

HD 4628 

Ks 

2010-08-03 

8.69 

HD 4676 

BrG 

2009-10-28 

7.9 

HD 4813 

BrG 

2010-08-05 

8.95 

HD 5448 

BrG 

2009-10-28 

8.89 

HD 7439 

BrG 

2010-08-05 

9.28 

HD 10307 

BrG 

2009-10-28 

8.23 

HD 10476 

BrG 

2010-08-05 

8.99 

HD 11171 

BrG 

2010-08-05 

8.21 

HD 11636 

BrG 

2010-08-05 

9.02 

HD 13161 

BrG 

2009-10-28 

9.2 

HD 13974 

BrG 

2009-10-28 

8.69 

HD 14055 

BrG 

2009-10-28 

8.72 

HD 16673 

BrG 

2010-08-05 

9.04 

HD 17093 

Ks 

2009-10-28 

8.61 


Lick IRCAL observations with no companions detected. The 
depth is the 5-<r limit reached at 4-5" range, where the contrast 
is highest. The full table is available online. 

we have identified 5 new companions. We discuss each new 
system below. 

HD 1404: This is an A2 star located 41 pc from Earth. 
We detect a companion 6.7 magnitudes fainter at Ks and 
Br -7 that over a three year period has remained at a pro¬ 
jected separation of 6 . 6 " with position angle (PA) of 144 
degrees. The companion has Ks~11.2, or Mk s ~8.1. At an 
age of 450 Myr (Vican 2012), the Baraffe et al. (1998) mod¬ 
els suggest a mass of ~0.17 Mq and effective temperature of 
~3200 K (spectral type M4) for this companion. This A+M 
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Figure 2. IRCAL 5 -a detection limits for the sample including 
detected known systems (green triangles) and new companions 
(red diamonds). In blue we show the median 5-<r sensitivity of 
our survey, along with the ±1 and ±2-<r ranges. 

binary is separated by 271 AU and hosts a debris disk of 
radius 22 AU around the primary A star (Thureau et al. 
2014). 

HD 168151: This system is an F star situated 32 pc 
from Earth. We detect a companion 3.3” from the primary 
in observations carried out in 2010 and 2012. The compan¬ 
ion is 6 magnitudes fainter than the primary suggesting an 
absolute Ma n magnitude of 7.4. The system’s age is esti¬ 
mated to be 1.3 Gyr (Vican 2012). At that age, the Baraffe 
et al. (1998) models predict a T~3400 K and 0.3 Mq (~M3) 
dwarf would have that absolute magnitude. 

HD 140538: This G-star had a previously known com¬ 
panion at a projected separation of 4.2" (Eggleton & 
Tokovinin 2008). With our Lick survey, we have resolved the 
companion into two equal brightness components separated 
by 0.22". At a distance of 15 pc, this separation corresponds 
to just over 3 AU. We estimate J, Fell, and Br -7 magnitudes 
of 8.5, 7.9, and 8.1 for the pair, with one component being 
~0.1 mag fainter than the other. The system has various 
age estimates, but as suggested by Vican (2012), we adopt 
the chromospheric age estimate of 3.6 Gyr. The Baraffe et al. 
(1998) models suggest a mass of 0.25-0.3 Mq (^M3) for each 
resolved component. 

HIP 4^^0: This is a late-K star 17 pc away that 
has been previously suggested to be an astrometric binary 
(Makarov & Kaplan 2005). We detect a companion 0.4" from 
the primary with Amag of 2.3-2.4 across J, H, and Ks. The 
JHKs magnitudes of the companion are 8.7, 8.2, and 7.9, 
respectively. The J-Ks color (0.8) is suggestive of an early 
to mid M-dwarf. No age estimate is available for HIP 42220 
(Vican 2012), but K-stars in our sample have median ages 
of 2 Gyr and dispersion of 2.9 Gyr. From the absolute Ks 
magnitude (6.7) and an age range of 100 Myr-5 Gyr, we es¬ 
timate the mass of the companion lies roughly in the range 
0.3-0.35Mq . The information listed in Makarov & Kaplan 
(2005) can be used with their Equation 2 to derive a lim¬ 
iting mass ratio for the unseen astrometric companion. In 
this case, assuming a circular orbit, total system mass of 
0.9-IMq , and minimum projection effects, we find our de- 


Figure 3. Disk detection frequency among the DEBRIS sample, 
divided by spectral type and binarity. 

tected companion is massive enough to be consistent with 
the astrometric signature observed by Makarov & Kaplan 
(2005). 

HD 110833: This system is a K3+K3 binary 15 pc away. 
Raghavan et al. (2010) list this as an equal mass spectro¬ 
scopic binary with a 271-day period. Vican (2012) estimates 
an age of 2.2 Gyr for this system. Our AO work reveals 
another star 1.6" from the primary with Amag of 4.2-4.4 
across J, H, and Ks. This projected separation amounts to 
~24 AU. The absolute JHKs magnitudes of the companion 
are approximately 8.7, 8.4, 8.1, and the Baraffe et al. (1998) 
models suggest effective temperatures of ~3200 K and mass 
~O.17M 0 for that age. The detected companion thus consti¬ 
tutes a mid-M dwarf in the system. 


4 DEBRIS DISKS AND MULTIPLICITY 

Systems in the DEBRIS sample have been observed with the 
Herschel spacecraft with PACS (100 and 160/rm; Poglitsch et 
al. 2010) and in some cases with SPIRE (250, 350, 500/rm; 
Griffin et al. 2010) and/or additionally at 70/160^nn with 
PACS. Furthermore, most also have observations with IRAS, 
Spitzer-MIPS, AKARI, and WISE. Hence, the spectral en¬ 
ergy distribution (SED) of these objects can be sampled in 
a broad wavelength range and any IR excess emission can 
be characterized and studied in detail. Full descriptions on 
the observation strategy, source extraction, SED modeling, 
and selected results can be found in Matthews et al. (2010) 
and other DEBRIS publications (e.g., Lestrade et al. 2012, 
Booth et al. 2013, Thureau et al. 2014, Matthews et al., in 
prep; and references therein). To summarize, stellar atmo¬ 
sphere models are first fit to data shortward of about 10 /rm 
to determine the photospheric flux. These predictions are 
then compared to the mid to far-IR photometry to look for 
the excesses indicative of emission from cool dust. Disk sys¬ 
tems are those with at least one photometric 3(J excess. To 
characterise the disk properties, a blackbody is fit to the 
star-subtracted SED and modified by a factor (Ao/A ) /3 for 
wavelengths longer than Aq to account for inefficient grain 
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emission at long wavelengths. For most cases, Ao is left as a 
free parameter, though in cases with limited photometry it is 
fixed (see details in Matthews et al, in prep). We note, how¬ 
ever, that a number of systems have poorly constrained tem¬ 
peratures or are required to have cold disks (Tdust <20 K), 
which could instead be a result of contamination from extra- 
galactic sources (Sibthorpe et al. 2013; Krivov et al. 2013; 
Caspar & Rieke 2014). For purposes of this paper, we in¬ 
clude these ~8-10 systems as part of our statistical analysis. 
Among the DEBRIS sample, we find 76 systems with IR ex¬ 
cess emission suggestive of circumstellar (or circumbinary) 
disks. In two cases (Fomalhaut and HD 223352), two stars in 
the same system each show IR excess which implies 78 disks 
among individual components the sample (see Section 4.4). 


4.1 Global Statistics 

Among the DEBRIS sample, there are 21 multiple systems 
which host circumstellar disks. There are several ways to 
approach the statistics of stellar multiplicity and disk inci¬ 
dence in our sample. We first consider the disk detection 
frequency among both single and multiple stars by spectral 
type. We note that our ability to detect a disk, or, more 
accurately, detect an infrared excess, will depend on the 
spectral type because our survey is flux-limited. The dis¬ 
tribution of distances for single and multiple star systems is 
similar so no distance biases are introduced. We note that 
multiple systems, regardless of the number of components, 
are counted only once. Considering the entirety of our 449- 
star sample, the 76 disk systems identified suggest an over¬ 
all detection frequency of 17 ± 2%. For single stars, this is 
55/261 = 21 ± 3% and for multiples, 21/188 = 11 ± 3%. 
Figure 3 summarizes the disk detection frequency in the 
DEBRIS sample as a function of spectral type and bina- 
rity. The detection of disks around A-stars is comparable 
regardless of multiplicity (35lg% for singles and 29tg% for 
multiples; see also Thureau et al. 2014). Among the lower 
mass stars, multiples host fewer disks than single star sys¬ 
tems. For the FGK sample, we find single stars have a disk 
detection frequency of 25lg% while multiple systems have 
8±1%. A two-sided Fisher test for the A-stars reveals a p- 
value of 0.64, which implies any different in detected disk 
frequency is not statistically significant. In contrast, for the 
FGK sample the Fisher test returns a p-value of 3.5 x 10 -4 
implying a statistically significant difference in the detected 
disk frequency between single and multiple star systems. 

An equivalent way of examining the frequency of disks 
around multiple-star systems is to consider the multiplic¬ 
ity fraction around dusty stars and those not bearing disks. 
This fraction is the ratio of the number of binary or multi¬ 
ple systems compared to the total number of systems in the 
disk or diskless samples; that is, the individual components 
within a multiple system are not considered. Figure 4 illus¬ 
trates the multiplicity fraction as a function of spectral type. 
Again, A-stars show similar multiplicity fractions regardless 
of whether we consider the disk and non-disk sample. Disk¬ 
bearing stars among the lower mass stars, however, have 
lower multiplicities. That is, it is rarer to find an FGK disk¬ 
bearing system that is also a binary or multiple. 

In addition to the incidence of infrared excesses among 
the DEBRIS sample, we examine properties of the disk, 
namely the distribution of dust temperature and the disk 



Spectral Type 


Figure 4. Multiplicity fraction among the DEBRIS sample, di¬ 
vided by spectral type and presence or absence of detected disks. 



Log Age (Gyr) 


Figure 5. Cumulative age distribution for A and FGK stars in 
both single and multiple star systems. 


luminosity divided by the stellar luminosity (Lm/Lbol', or 
fractional luminosity), to see if any differences or trends are 
evident in the sample. A KS test reveals no difference be¬ 
tween single and multiple stars in terms of dust temperature 
(p=0.9) or Lm/Lboi (p=0.5), as was also demonstrated for 
A stars by Thureau et al. (2014). This holds regardless of 
whether we consider the full sample or break it as A or 
FGK stars: for dust temperatures we find p=0.3 for A sin¬ 
gles vs multiples and similarly p=0.9 for FGK systems. For 
Lm/Lboi, we find p=0.3 and 0.4 for A and FGK systems, 
respectively. For systems with detected disks, the basic prop¬ 
erties of these disks don’t appear to be correlated with the 
multiplicity of the system. 


4.2 Disk Frequency Around A and FGK Stars 

The frequency of disks around A stars, as well as their prop¬ 
erties, is the same regardless of multiplicity. The situation is 
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different for FGK stars in that while the disk properties are 
similar, the frequency of disks is lower among the multiple 
systems. 

The dynamical effects of a second star can serve to dis¬ 
rupt the disk, accelerating its dispersal. As such, the disk de¬ 
tection frequency among multiples should be lower, or equiv¬ 
alently, disk-bearing systems should more likely be single 
stars. The discrepancy between the A-star and FGK sam¬ 
ple, however, could be a result of a selection bias due to their 
ages. Examining ages from Vican (2012), we find dusty A- 
stars in the DEBRIS sample have a median age of about 
~0.3 Gyr, whereas the dusty FGK-stars are older, at a me¬ 
dian age of ~4 Gyr. On the other hand, the binary fraction 
of disk-bearing systems for <1 Gyr and >1 Gyr FGK stars 
is 29lg 4 % and 21+g%, respectively. While we expect older 
FGK stars to be less likely to host disks, the multiplicity 
difference we observe is not statistically significant. 

We also considered the age distribution of A, F, G, and 
K stars for both single and multiple stars regardless of the 
presence of disks. If the ages were different as a function of 
multiplicity, these could account for the observed difference 
in the disk frequency. However, both single and multiple 
star distributions look very similar (see Figure 5). The KS 
test does not show evidence that singles and binaries of any 
spectral type are drawn from separate age distributions. 

As a further test to ascertain the difference between 
A and FGK binaries with disks, we constructed samples of 
systems eliminating either very wide binaries (>500 AU) or 
very close binaries (<1 AU). In both cases, these compan¬ 
ions tend to be more common for A stars than lower-mass 
stars and could contribute to the discrepancy (see Figure 2 
in Duchene & Kraus 2013). These systems have little effect 
on the disks Herschel is sensitive to. However, in our sample 
there are only a handful of these wide or close systems and as 
such the statistics do not appreciably change when eliminat¬ 
ing these. In fact, the distribution of periods or semi-major 
axes between A and FGK multiples among this sample is 
not appreciably different. 

The discrepancy between A and FGK disk-bearing bi¬ 
naries does not appear to be a result of binary properties 
or an age bias in our sample. A possible explanation is that 
we are missing a significant number of companions among 
A-stars, as suggested by the higher multiplicity fraction 
(>50%) in Duchene & Kraus (2013). However, the A-stars 
in the DEBRIS sample are generally well-studied and we 
would require 6-7 additional binaries, all without detected 
disks, to yield a comparable difference between disk-bearing 
and non-disk bearing systems as observed among the FGK 
stars. 

4.3 Disks in Binary and Multiple Systems 

As previously mentioned, the binary period distribution of 
the full DEBRIS sample is consistent with prior studies of 
stellar multiplicity (Eggleton & Tokovinin 2008; Duquennoy 
& Mayor 1991; Raghavan et al. 2010). Figure 1 highlights 
the period distribution of disk-bearing systems, which ap¬ 
pears flat in comparison to the full sample. Figure 6 shows 
the cumulative distribution of periods for the non-disk and 
disk binaries compared to the distribution of Raghavan et 
al. (2010) and a flat distribution across all periods. At a 
glance, among the relatively low number of disk-bearing bi- 



0.0l— i -' ..i - i - .-i -'-■- 

-2 0 2 4 6 8 10 

Log Period (d) 


Figure 6. Cumulative distribution of calculated and measured 
orbital periods for multiples. Black shows the full sample, while 
blue are disk-bearing binaries. The Raghavan et al. (2010) distri¬ 
bution is shown as a dashed line and a flat distribution is shown 
with the dotted line. 


naries there appear to be fewer of these with periods between 
10 3 and 10° days. However, the period distribution of disk¬ 
bearing DEBRIS systems does not differ significantly from 
that of the non-disk sample. A KS test returns a p-value 
of 0.09 when comparing these two distributions. For us to 
assert that the two samples to are drawn from different dis¬ 
tributions we require a p-value of 0.05 or lower. Hence, we 
cannot claim that the distribution of disk and non-disk bina¬ 
ries differ from each other or from the general distribution of 
Raghavan et al. (2010) or from a flat, uniform distribution. 
While the various distributions look tantalizingly different 
at first glance, we cannot rule out they are drawn from the 
same underlying distribution with the limited statistics of¬ 
fered by this survey. 

Another way to examine the disk-bearing binaries is by 
examining the stellar separation. Figure 7 shows the stel¬ 
lar separation compared to the location of the dust in the 
system for multiples in this work, as well as others from Ro¬ 
driguez & Zuckerman (2012) and Trilling et al. (2007). For 
the dust, we plot the semimajor axis of the dust assuming it 
is composed of large, blackbody grains. Some disks are spa¬ 
tially resolved and discussed elsewhere (e.g., Kennedy et al. 
2012a,b; Booth et al. 2013). For A-stars, the disk location for 
resolved systems tends to be a factor of ~l-2.5 times that es¬ 
timated by the assumption of blackbody grains; FGK disks 
may be larger (see Pawellek et al. (2014)). As such, some of 
the systems may be shifted rightward in the plot. For clarity 
and consistency, we use only the dust semimajor axis as de¬ 
rived from the SED rather than any resolved radius. Some 
systems may be shifted upward because of projection effects 
(or downward if observed near periastron) in the plot, but 
this tends to be a small correction (e.g., Dupuy & Liu 2011). 
Some multiples are connected by dotted lines corresponding 
to the stellar separations. For example, a triple may con¬ 
sist of a close binary surrounded by a dust disk and a more 
distant stellar companion. Both the separation of the close 
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Figure 7. Stellar separation vs blackbody dust semimajor axes 
for the dusty multiple systems in the DEBRIS sample and from 
other surveys in the literature. The grey denotes the approximate 
region at which gravitational effects can perturb the disk. Triples 
where the dust lies between the AB and C pair are connected 
with a dotted line. HD 223352 is the system with the slanted line 
as the system hosts two separate debris disks. With stellar sepa¬ 
rations exceeding 150 kAU, the three components of HD 216956 
(Fomalhaut) are beyond the plot range in this figure. 


binary and the more distant companion would be plotted 
and connected in Figure 7. 

Figure 7 also shows a grey region representing the area 
where the gravitational influence of a companion would dis¬ 
rupt the disk. This is done by computing the critical semima¬ 
jor axis using the relationships in Holman & Wiegert (1999). 
This is the distance at which a test particle would survive 
for less than 10 4 times the binary orbital period. For these 
relationships, we adopt 0.5 for the mass ratio and 0.4 for the 
eccentricity and derive critical semimajor axes values of 0.15 
and 3.4 times the stellar semimajor axis (see Tables 3 and 
7 in Holman & Wiegert 1999 for values for other mass ra¬ 
tios and eccentricities). Systems located in this region would 
have their disks quickly cleared out by the stellar compan¬ 
ion. Of the 26 DEBRIS components plotted, 4 (HD 223352, 
99 Her, HIP 14954, HIP 73695) lie in this unstable region, or 
lb±l°%. If we generate a random sample populating the di¬ 
agram, either with a uniform distribution or a distribution 
following the period distribution of Duquennoy & Mayor 
(1991) and with a disk population uniformly spread between 
disk radii of 0.12 to 1000 AU, we find that ~20% of compo¬ 
nents lie in the unstable area. Errors due to unknown incli¬ 
nations or dust grains deviating from blackbodies were not 
included in this simulation, but their effects are expected to 
be minor (see above and Dupuy & Liu 2011). This suggests 
the location of the dust in the system may not be strongly 
correlated with the location of the stellar companion. 

As previously mentioned, the fractional luminosity, 
Lir/L boi , of multiples with detected disks is not significantly 
different from that of single stars within the DEBRIS sam¬ 
ple. Figure 8 compares this fractional luminosity against the 
ratio of the dust semimajor axis and the stellar separation. 
Systems located towards the left of the plot are circumstel- 


Figure 8. Ratio of dust to stellar separation compared with 
Lir/L* for dusty multiple systems. The grey denotes the ap¬ 
proximate region at which gravitational effects can perturb the 
disk. As before, triples are connected and HD 223352 is the sys¬ 
tem connected with a slanted line as two separate debris disks 
are present in the system. HD 216956 (Fomalhaut) has a ratio 
< 10 —3 and is not shown in this figure. The system with frac¬ 
tional luminosity ~4% is BD+20 307. 

lar in nature where the dust orbits a single star, yet there 
is a distant stellar companion in the system. Systems to¬ 
wards the right are circumbinary in nature in that the dust 
surrounds a pair of stars. Again, the unstable zone is high¬ 
lighted as discussed above. The spread in fractional lumi¬ 
nosity is comparable between the circumstellar disks (dust 
location / stellar separation ratio < 0.1) and the circumbi¬ 
nary disks (ratio > 3). Combined with the results of Fig¬ 
ures 6 and 7, this suggests the properties of disks around 
disk-bearing binary stars do not, in general, depend strongly 
on the orbital properties of the binary. 

4.4 Individual Binary Debris-Disk Systems 

In this section, we highlight a handful of disks around binary 
or multiple star systems. These stand out in our sample as 
we describe below. 

HD 223352: As detailed in Phillips (2011), HD 223352 
is a quadruple system in which two components, A and C 
(HD 223340), host circumstellar disks. The AB pair is sepa¬ 
rated by 3.9” (164 AU) and the C component is much farther 
away at 74” (3100 AU). Recent observations have resolved 
the B component as two separate stars (De Rosa et al. 2011). 
The two disks in the system are located around the A and 
C stars, with estimated blackbody dust locations of 27 and 
10 AU. The ratio of the dust to stellar separation for the 
AB pair is 27/164 = 0.16 placing it just inside the unstable 
region in Figures 7 and 8, whose upper bound is set at a 
ratio of 0.15. For Figures 7 and 8, we have taken care to 
connect the two disk-bearing components together as each 
disk is associated with a different star separation within the 
system. 

Very few debris disk systems are known in which two 
stars host individual disks. Disregarding ~10-20 Myr-old 
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(and younger) systems whose disks may be primordial in na¬ 
ture (see, e.g., Rodriguez et al. 2014, and references therein), 
the only other multiple system known to host two disks is 
that of Fomalhaut, also part of this survey (see Kennedy et 
al. 2014). As a member of AB Dor (Malo et al. 2013; Zuck- 
erman et al. 2011), HD 223352 has an age of ~100 Myr and 
represents a rare opportunity to study the development and 
evolution of disks around a pair of co-eval stars. 

HD 165908: Commonly known as 99 Her, this system 
has been described in detail in Kennedy et al. (2012a). With 
a stellar separation of ~16 AU and a blackbody disk dust 
location of ~47 AU, it is one of the unstable disk systems in 
our study. However, the disk has been resolved by Herschel 
to be located ~120 AU from the binary pair, which places 
the system beyond the grey unstable zone in Figure 7. A 
more detailed discussion of the disk and its interaction with 
the binary can be found in Kennedy et al. (2012a). 

HIP 1^95^: An F8 primary orbited by a widely sepa¬ 
rated M star at ~223 AU. Like 99 Her above, HIP 14954 
appears to have dust located in an unstable configuration 
given the blackbody dust location of ~164 AU is compara¬ 
ble to the stellar separation. Roberts et al. (2011) provides 
orbital parameters for this system: a period of 2029 years, 
a semi-major axis of 9.87” (223 AU), and eccentricity of 
0.26. Using these orbital elements, the relationships of Hol¬ 
man & Wiegert (1999) predict critical semimajor axes of 59 
and 708 AU. That is, orbits whose semimajor axis is smaller 
than 59 AU (or larger than 708 AU) are stable on timescales 
longer than 10 4 times the stellar period (~20 million years 
in this case). The fact that the dust is located beyond 59 AU 
suggests it should have been quickly disrupted or instead has 
been recently produced. However, this system shows possi¬ 
ble confusion in the Herschel data and the infrared excess 
could be originating from a background source. 

HIP 73695: A triple system in which the primary is or¬ 
bited by a 0.3-day binary. The binary orbits the primary 
with a period of 206 years and eccentricity 0.55 (Eggleton & 
Tokovinin 2008); this amounts to a semi-major axis of about 
37 AU. The dust is around the primary at a blackbody dust 
location of ~11 AU. The location of the dust is uncertain 
as the infrared excess is detected only at 70/rm. With the 
binary properties listed in Eggleton & Tokovinin (2008), we 
estimate critical semimajor axes of 4 and 138 AU. As with 
HIP 14954, the dust appears to be located beyond the in¬ 
ner critical semimajor axis. Additional information would 
be needed to make a definitive statement on the stability of 
the disk in the system. 

5 CONCLUSIONS 

The DEBRIS sample consists of 449 AFGKM systems ob¬ 
served with Herschel to search for circumstellar disks. In 
this study, we have examined these stars for signatures of 
stellar multiplicity. In this sample, 42% of the targets are 
binary and higher order stellar multiples, as noted by others 
in the literature and our own adaptive optics observations 
(Table 3). This sample allows us to examine the influence of 
stellar multiplicity on the incidence and properties of debris 
disks. 

Disk systems among the sample are summarized in 
Matthews et al. (in prep), Thureau et al. (2014), and 


Sibthorpe et al. (in prep). Despite the large number of sys¬ 
tems considered, the number of disk-bearing systems (78 
across all spectral types) is lower than considered in Ro¬ 
driguez & Zuckerman (2012). This is a result of the way 
these two samples were constructed. Nevertheless, an ad¬ 
vantage to this study is that it allows us to consider both 
the disk frequency (as in Trilling et al. 2007) and the bi¬ 
nary/multiple frequency (as in Rodriguez & Zuckerman 
2012) simultaneously for the same sample. Our results are, 
at a glance, similar to both prior studies. However, an ex¬ 
amination of the period distribution of the multiples in the 
sample shows no statistically significant difference with re¬ 
gards to whether or not these systems host debris disks and 
instead follows more closely the distribution previously re¬ 
ported for stellar multiples in general (e.g., Duquennoy & 
Mayor 1991; Raghavan et al. 2010). 

Unlike prior studies, the larger sample allows us to 
break up the results as a function of spectral type. Among 
A-type stars, as primarily examined in Trilling et al. (2007), 
Rodriguez & Zuckerman (2012), and Thureau et al. (2014), 
disks are detected just as frequently among single stars as 
in multiple star systems. Equivalently, disk and non-disk A- 
star systems show comparable multiplicity. The results for 
FGK stars, however, differ substantially from the A stars. 
In FGK stars, the presence of a second star reduces the 
likelihood to find a circumstellar or circumbinary disk in 
the system. We have ruled out potential biases due to age 
or binary properties, but have been unable to find a con¬ 
vincing mechanism for this difference compared to A stars. 
Having a significant number of missing binaries among A- 
stars remains a plausible explanation. A larger sample of 
well-characterized FGK stars may be needed to further ex¬ 
amine the disk detection frequency as a function of stellar 
multiplicity. 

Among the disk-bearing multiples, both circumstellar 
and circumbinary disk systems exist. The basic properties 
of detected disks in binary systems, namely the tempera¬ 
ture of the dust and the fractional luminosity, are compa¬ 
rable to those in single star systems. There are a few sys¬ 
tems (described in more detail in Section 4.4) in which the 
stellar separation and the location of the dust (assuming 
large blackbody grains) are comparable to each other. This 
would result in an unstable scenario in which the gravity of 
the stellar companion would readily disrupt the disk. How¬ 
ever, limitations in our knowledge of the orbital parameters, 
namely the projection on the sky, as well as the possibility 
of inefficient grain emission in the disk, can conspire to make 
a system appear to have an unstable configuration when it 
does not. 

We interpret these results as follows. Any sample of bi¬ 
nary stars dominated by FGK stars (as is this one) will in 
general demonstrate the period distribution of stellar mul¬ 
tiples of Duquennoy & Mayor (1991) and Raghavan et al. 
(2010). The distribution of such a sample will not vary sig¬ 
nificantly between systems with detected dust disks and the 
general population. A binary star hosting a circumstellar 
or circumbinary disk will exhibit properties consistent with 
other (non-disk bearing) binaries and the dust properties 
in the system will be comparable to those of single stars 
with disks. However, the gravitational influence of the com¬ 
panion star(s) will accelerate the collisional evolution of the 
disk. As such, a sample of binary stars will have a detectable 
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excess for a shorter period of time. Hence, disks will be less 
readily detected among a sample of old binary or multi¬ 
ple stars than single stars of comparable spectral types and 
ages. However, given the stochastic nature of debris disk 
formation as a function of time, this signature is diminished 
among samples dominated by stars with ages larger than a 
few hundred million years. 

The increasing number of planets found orbiting pairs 
of stars shows that planet formation in binary systems is 
not an uncommon product of planetary evolution. Binary 
and multiple star systems are plentiful and our work sug¬ 
gests that the processes behind disk and planet formation 
are similar between single and multiple star systems. 
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